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Qh Abstract 



Fractional integral operators connected with real- valued scalar functions of matrix argument are applied 
in problems of mathematics, statistics and natural sciences. In this article we start considering the case 
of a Gauss hypergeometric function with the argument being a rectangular matrix. Subsequently some 
fractional integral operators are introduced which complement these results available on fractional operators 
in the matrix variate cases. Several properties and limiting forms are derived. Then the pathway idea is 
incorporated to move among several different functional forms. When these are used as models for problems 
in the natural sciences then these can cover the ideal situations, neighborhoods, in between stages and paths 
leading to optimal situations. 
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1. Introduction 



The importance of anomalous reaction/relaxation and transport/diffusion is well recognized in many 
disciplines including physics, chemistry, biology, and engineering. Despite this fact, anomalous relaxation 
and transport are not well understood, and there is the need to develop mathematical and statistical models 
with predictive power. Recent developments of fractional calculus in terms of integro-diffcrcntial operators 
that provide a unifying framework to model key aspects of anomalous relaxation and transport like non- 
locality, non-Markovian (memory) effects, and non-Gaussian (Levy) processes became available (cf. Mathai 
et al. [11], Nair [4]). The extension of such developments to matrix-variate statistical densities in general 
(Mathai and Haubold [8], Mainardi [13]) and to the specifically interesting case of Mittag-Leffler functions 
and matrix-variate analogues (Mathai [7]) have been achieved more recently. Such results are applicable 
to the solution of linear coupled fractional differential equations (Lim et al. [1]) as well as to the handling 
of fractional Poisson probability distributions (Laskin [12]). Further use of fractional integral operators 
connected with real- valued scalar functions of matrix argument can be utilized for fractional matrix calculus 
(Phillips [3]), numerical solution of fractional diffusion- wave equations (Garg and Manohar [6]), stability 
analysis of fractional-order systems (Jiao and Chen [9]), and probably for the application of the matrix- 
variate Mellin transform in radar image processing ( Anfinsen and Eltoft [10]). 

Let X — (xij) be a p x r, r > p matrix of real distinct scalar variables x^-'s. Let A be a p x p real 
positive definite constant matrix, that is, A = A' > O, prime denoting the transpose. Let B be a constant 
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positive definite r x r matrix, B > 0. Let A? and denote the positive definite square roots of A and 
B respectively. Let X be of full rank p. Then Z x = A^XBX' 'A* is symmetric positive definite matrix. 
In this paper we will consider only real matrices. The corresponding results in the complex domain can be 
done parallel to the real case. The following standard notations will be used. Real matrix variate gamma 
will be denoted by 

T p (a) = n P -^T(a)T(a l)..T(a - J -=^), (1.1) 

for 5R(a) > where 3?(-) denotes the real part of (•). It can be shown that T p (a) has the integral 
representation 

r p (a) = / (Sp-^e-^dS, K(a) > ^— - (1.2) 



J5>0 ^ 

where |5| means the determinant of the p x p positive definite matrix S and tr(S) denotes the trace of S. 
The wedge product of the p( - P 2 1 " 1 ^ differentials dxij 's will be denoted by 

dS = JjAdsy, A= wedge (1.4) 

i>j 

The type-1 beta integral is given by 

« = TTT^ = / ISI-^lz-sr^dS, (1.5) 

ip(a + pj Jo<x<i 

for 3?(a) > 2^-, 5ft(/3) > 2=^, where 5 is p x p positive definite and O < S < I means S>0,I-S>Oot 
all the eigenvalues of S are in the open interval (0,1). In general, J x means the integral over X. The type-2 
beta integral is given by 



B p (a, P) = [ ISl^^lI + S\-( a+ ^dS 
Js>o 
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\U\ fi -*P \I + U\- (a+fj Uu, »(a) > ^— 5R(/3) > - (1.6) 



/(7>0 

Let 

Z x = AJIBX'AJ and Zy = A^By'A^ 

where X = (xij) and Y = (yij) are p x r,r > p matrices of real elements, and of full rank p. Consider the 
evaluation of the integral 

h= [ \Z x \ a \I-Z x \ c - a - P -¥\I-Z Y Z x \- b dX, (1.7) 

for O < Zx < 1,0 < Zy < I. This integral corresponds to the Euler integral for Gauss hypergeometric 
function. Let 

U = A^XB^ => dU = \A\%\B\%dX (see Mathai, 1997, equation (1.2.20)) (1.8) 
Then, after integration over the Stiefel manifold, 

rp 

Z X = UU' = V ^dU =JL±-\v\i- t S L dV (1.9) 
(see Mathai, 1997, Theorem 2.16 and Remark 2.13). Now, the integral I\ becomes, 



7T 2 



r P ( 



P\ If 

+ i p+i 



x f \v\ a+ i-^\I- V\ c - a -^\I - Z Y VZ Y \- b dV. (1.10) 
Jo<v<i 



Let us expand the last factor in the integrand in terms of zonal polynomials. For a discussion of zonal 
polynomials see Mathai et al. (1995). 

\I Z Y VZ Y \- b = f] £ ^C K (Z Y V) (1.11) 

fc=0 AT 

where Ck(-) is the zonal polynomial of order k, K — (ki,k 2} fc p ), fci + ... + k p — k and 



- H(6 ^— 



r p (&,#) = J] r(6 + % - = r p (6)(6)x (1.12) 



i=i 

and (a) m is the Pochhammer symbol 

(a) m = a(a + I). ..(a + m - 1), (a) = l,a^0. 

We can evaluate the integral 

f wrt-^ii vr^c K (z Y v)dv = Tp{a y?'*?i c - a) c K (z Y ) a, is) 

Jo<v<i r p (c+2,-K) 

(see Mathai, 1997, (5.1.26)). Note that interchange of integrals and sums is valid here. Then 

11111 r p (§) r p ( c + |) 
^EE^f^. (i.u, 

2. Hypergeometric Functions of Rectangular Matrix Argument 

In the notations of hypergeometric function of matrix argument the series part of (1.14) can be written 
as a 2-Fi- That is, 

y> \ - OkO + \)k Ck(Z y ) , r r 

2^1^ fr+ i W ^—=2F 1 (a+-,b;c+-;Z Y ), \\Z Y \\ < 1 (2.1) 

where ||(-)|| denotes a norm of (•). Hence we have the following theorem: 
Theorem 2.1. For X,Y 7 A,B as defined in Section 1 

r.„, |A|S|s|$r p (§)r p (c + §) 



2 F 1 (a+-,b;c+-;Z Y ) 



2 2' 7r 2 r p (a + §)r p (c - a) 

x / \Z x \ a \I - Zx^-^ll - Z Y Z x \- b dX (2.2) 

for O < Z x < 1,0 < Z Y < J,R(c-a) > ^.^(a) > -§ + ^ 

For establishing some limiting forms and pathways we need the following results, which will be stated 
as lemmas. 
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Lemma 2.1. For K = (fci, k p ), ki + ... + k p = k, (a) K = H^ii - - L T')k j , 



lim( g -l) fe (-i-)* = l. 
g->l q — 1 



Note that 



Lemma 2.2. 



lim |J+(t?- =c - tr (^). (2.4) 

ejr->l 



Proof: Let Xj, j = 1, ...,p be the eigenvalues of Zy. Then 

p 

|7+( g -l)Z K | = J](l + ( g -l)A,). 

But 



n[lim{l + (g - 1)Aj}-^t] = ]J e"^ 1 = e"^=i Aj 

= e- tr(Zy) . (2.5) 



3. Fractional Integral Operators 

Let and Z Y be as defined in Section 1. Let the generalized fractional integral operator of matrix 
argument be defined and denoted as 

( D x a .f)(X) = j \Z X - Z Y \ a - E i 1 f(Z Y )dY 

L p( a ) Jz x >z Y >o 

\Z X \ a - V -^ f .„_£+! 



L» 



/ \I-Z x i Z Y Z x i r E i 1 f(Z Y )dY. 

JZ x >Zy>0 



Make the transformations U = A2YB2 , V = UU', W = Z X 2 VZ X 2 . Then integrating out over the Stiefel 
manifold we have 



oD x a f 



\A\t\B\iT p (a)T p (%) 



x f \I -W\ a -^ f{Z x WZ x )\W\^^AW. (3.1) 

Jo<w<i 



Consider the special cases of /(•). Consider the operator operating on a power function. 
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Case 1. Let f(Z x ) = \Z X \ V - Then 



o^x \ A x\ 



Jo<w<i 



\A\t\B\*r p (a)T p (%) Jo<w<i 

\Z x \^+n-^^ T^+n) ^r +ri)> p-l (32) 



L4|i|5|5r p (§) ^(a+i+r/)' ^2 



Case 2. Zonal polynomial of order k, f(Zx) — Ck(Zx) 
Then going through the same steps as above 



ry | Q +Il_£±i IE 

ZX 2 2 7T 2 I r p+1 



\A\2\B\2T p (a)T p (^) Jo<w<i 
x \I - W^' 2 ^ Ck(Z x W Z x )dW, (Mathai, 1997, equation (5.1.26)) 



\A\i\B\ P 2T p {a)T p ^) T p (a + ^,K) 

\Zx\^2ZJ^lZ (i)KC K (Z x ) 
|A|5|B|5r p (a+§) (<*+§)*• 



4. Extended Saigo Operators 

Let 

= 2 Fi(a,6;c;/-Z x ^Z y Z x i ) ( / ) (Z y ),Z x > Zy. 

Here 2 F\ is a Gauss hypergeometric function of matrix argument I — Z x 2 ZyZ x 2 and (j)(Zy) is an arbitrary 
function so that 

/ \I-Z x i Z Y Z x i \ a - £ ^C K (I-Z x i Z Y Z x i ) ( f>{Z Y )dY <oo. 

Opening up the 2-F1 in terms of zonal polynomials and then substituting U = A^YBA? , V" = £/£/' we have 



oD^^a, b-c-I- Z x *ZyZ x *)4>(Z y ) 
(a) K (b) K \Z x \ a - P -^^ 



EE- 
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k=0 K kl(c) K \A\i\B\iT p (a)T p ^) 

i - z^vz x ^\ a - v -^c K {i - z x h 'Z x ^)<i>{V)\V\i- v -^ dV. (4.1) 



o<v<z x 



Consider the special case <fi(V) = \V\ n and then make the transformation Z X 2 VZ X 2 — W. Then the right 
side of (4.1) reduces to the following: 

£^ k\(c) K \A\i\B\iT p {a)T p ^)Jo<w<i 
The integral part reduces to the following: 
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f iwr+i- 1 ^ 1 ]! -W^-^Cnil -W)dW 

Jo<w<i 

= [ |T| Q - 2 * i |/-T|" + i- JE * 1 C , A-(T)(iT 

JO<T<I 



r p (a,iQr p (T? + §; 

r p (a)r p (?7+ |) ( a )jf 



Or (J) 



r p (a + ?? + §) (a + ry+l)^ 
Substituting back and denoting the left side by Jx we have 



X \A\i\B\ir p (Z)T p (a + r,+ %y X] 

r 

x 3 F 2 (a,b 7 a;c,a + f]+ -;/). 



5. Some Statistical Considerations 

i i i i 

Let U\ = AlXiBiX[Al , [7 2 = A^XiBiX' 2 A 2 . Let X\ and X 2 be independently and exponentially 

distributed random matrices where A\ and A 2 are p x p real symmetric positive definite constant matrices, 

Bi is ri x n, ri > p and B 2 is ^2 x r 2, r 2 > P constant positive definite matrices, Xi is p x ri and X 2 is 

p x r 2 matrices of distinct real scalar random variables and the matrices be of full rank p, where a prime 

denotes the transpose and the square roots are the unique positive definite square roots of A x and A 2 

respectively. Let the densities of X\ and X 2 be denoted by fi{Xi) and f 2 (X 2 ) and the joint density, denoted 

by f(Xi,X 2 ) = fi(Xi)f 2 (X 2 ) due to independence. Further, let 

f j (X j )dX j = eye tr(A i x i"J x J A j\lX J .j = 1, 2 (5.1) 



where Cj is the normalizing constant. Make the transformations Vj = A? XjB? => dXj = \ Aj\ 2 \Bj\ ^dVj 



where dX,- is the wedge product of all differentials in Xk = ( x ij^) or <LYfc = nf=i Fljii Ada:y\ A = wedge 
product. Let Wj — VjV- and integrate out over the Stiefel manifold. Then we have 

dV j = ^-AW j \^- E i 1 dW j (5.2) 
1 pi 2 > 

and tr(-) denotes the trace of (•). Hence the rectangular matrix Xj having the exponential density means 
Wj is having a real matrix- variate gamma density and further, 

dX 3 = 2 dW . (5 .3) 

\Aj\+\Bj\iT p ft) 

If the densities of Vj and Wj are denoted by gj(Vj) and hj(Wj) respectively then 



hj(Wj)dWj = ,^^\Wj\^-^e-^dWj, Wj > O (5.4) 
1 p\ 2 > 



1 ,H_£±i 

nMWi = 4r^" <l,,v ' ,dl 'j < 5 - 5 > 

./•«.V,.d.V, = ^tjf* ^ ^ A h^ A h d Xj. (5.6) 

7T 2 
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5.1. Density of the sum 

Let us examine the density of the sum U\ + U\. The joint density of X\ and X 2 , denoted by f(X\, X 2 ), 
is given by 

f(X 1 ,X 2 )dX 1 AdX 2 = f 1 (X 1 )f 2 (X 2 )dX 1 AdX 2 . 
Make the transformations Vj = A]XjBj,Wj = VjVj then 



f(X 1 ,X 2 )dX 1 AdX 2 = /i(W)/ 2 (F 2 y 2 '){n \Aj\-Z\Bj\SydViAdV2 

3=1 

2 li 

= h(W 1 )f 2 (W 2 ){]J\A 3 \-^\B J \-i-^\W^-^}dW 1 AdW 2 . 
,-1 l p\ 2 ) 



Then U\ + U 2 = W\ + W 2 . Make the transformation U — W\ + W 2 and V = W\, the Jacobian is unity and 
Wi — V, W 2 = U — V and the integration is over the positive definite matrices U and V with U — V > O. If 
Xi and X 2 are independently distributed as in (5.1) then denoting the marginal density of U by f*(U) we 
have 

f*(U)dU= f fi(U-V)f 2 (V)dV (5.7) 
Ju>v>o 



= FT^ I^ - V^-^lV^-^e-^-^-^dV A dV. (5. 
Make the transformation Z = U~^VU~^ for fixed U, then 

/*(t/)dt/= " / |Z|^-^|/-Z|^-^dZAdC/ 

-d(7. 



r p (^) 

Thus the sum U is real matrix-variate gamma distributed. Note that the integral in (5.7) is Riemann- 
Liouville left-sided fractional integral for the power function or where the arbitrary function is of the form 
/(T) = |T|*-4 i . 
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